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ABSTRACT 

Context. The disk-outflow connection is thought to play a key role in extracting excess angular momentum from a forming proto-star. 
Though jet rotation has been observed in a few objects, no rotation of molecular outflows has been unambiguously reported so far. 
Aims. We report new millimeter-interferometric observations of the edge-on T Tauri star - disk system in the isolated Bok globule 
CB 26. The aim of these observations was to study the disk-outflow relation in this 1 Myr old low-mass young stellar object. 
Methods. The IRAM PdBI array was used to observe ^^C0(2-l) at 1.3 mm in two configurations, resulting in spectral line maps with 
1.5'' resolution. We use an empirical parameterized steady-state outflow model combined with 2-D line radiative transfer calculations 
and -minimization in parameter space to derive a best-fit model and constrain parameters of the outflow. 

Results. The data reveal a previously undiscovered collimated bipolar molecular outflow of total length ^ 2000 AU, escaping perpen- 
dicular to the plane of the disk. We find peculiar kinematic signatures that suggest the outflow is rotating with the same orientation as 
the disk. However, we could not ultimately exclude jet precession or two misaligned flows as possible origin of the observed peculiar 
velocity field. There is indirect indication that the embedded driving source is a binary system, which, together with the youth of the 
source, could provide the clue to the observed kinematic features of the outflow. 

Conclusions. CB 26 is so far the most promising source to study the rotation of a molecular outflow. Assuming that the outflow is 
rotating, we compute and compare masses, mass flux, angular momenta, and angular momentum flux of disk and outflow and derive 
disk dispersal timescales of 0.5 ... 1 Myr, comparable to the age of the system. 

Key words, circumstellar matter — ISM: jets and outflows — Stars: pre-main sequence — planetary systems: protoplanetary disks 
— ISM: molecules 



1. Introduction 

■ A comprehensive understanding of the complex relations be- 
tween a collapsing molecular cloud core, the forming young 
(proto-)star, its circumstellar disk, and the expelled bipolar jets 

[ is a necessary prerequisite to understand the formation of stars 
and planetary systems. In particular, the initial amount, evolu- 
tion, and re-distribution of angular momentum appear to be key 
factors in many of the processes which determine the protostel- 
lar fragmentation, the final stellar mass(es), and the existence 
and morphology of planetary systems. 

Rotation of protostellar jets and molecular outflows about 
their symmetry axis is believed to be a key mechanism to re- 
move angular momentum from forming protostars and, thus, 
at least partially, to solve the angu lar momentum problem in 
star formation (e.g.. iRay et al.ll2QQ7l) . As a consequence, many 
groups have searched for rotation signatures in jets and outflows. 



Send offprint requests to: Ralf Launhardt, e-mail: rl@mpia.de 

* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). Also based on observations collected at 
the Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, oper- 
ated jointly by the Max-Planck Institut fur Astronomic and the Instituto 
de Astrofisica de Andalucia (CSIC). 



Indeed, during the past few years, radial velocit}Q gradients 
have been observed across a number of protostell ar jets using 
high- resolution HS T spectra (IBac ciotti et al. 2QQ3: ICoffey et al.l 
I2QQ4I: IWoitas et al.ll2005l: fCoffev et a l. 2007). These were inter- 
preted as jet rotation tracing the helicoidal structure of the mag- 
netic field in accretion-ejection structures. As these observations 
are very difficult and the rotation signatures are weak, other 
groups have also searched for rotation signatures in molecular 
outflows. Whi le IPetv e t al. ( 2006) did not detect outflow rota- 
tion in HH 30, iLee et a l. ( 2007) report radial velocity gradients 
across SiO knots in HH 211 which they interpret as jet rotation. 
However, recent PdBI observations could not confirm this claim 
(Gueth et al. in prep.). 

In order to study the chemical composition, kinematics, and 
evolution of protoplanetary disks, we started the CID project 
("Chemistry In Disks"), a molecular line survey with the IRAM 
Plateau de Bure Interferometer (PdBI) of well-known disks 
such as those surr ounding LkCal5, MWC480, or DMTau 
dPutrev etal.1 l2QQ7h . Our main goal is to map, at high angu- 



^ Throughout this paper, we use the term "radial velocity" when we 
refer to the line-of- sight velocity in the coordinate system of the ob- 
server. For the radial component of the gas motion in the coordinate 
system of the source (i.e., gas moving radially away from the source), 
we use the term "radial outflow velocity". 



2 



Launhardt et al.,: Rotating molecular outflow in CB 26 



lar resolution, various molecular lines with enough sensitivity 
to derive molecular abundance variations versus radius. For the 
youngest source in our sample, CB 26, we also mapped the 
molecular outflow with the goal to study the disk -outflow re- 
lation. 

CB 26 (L 1439) is a small cometary- shaped Bok globule lo- 
cated 10° north of the Taurus- Auriga dark cloud, at a distance 
of ^ 140pc. A dense core with signatures of star formation is lo- 
cated at the south-west rim of the globule (iLaunhardt & HenningI 
OVRO observations of the mm dust continuum emis- 
sion and of the CO (1-0) line have revealed a nearly edge- 
on c ircumstellar disk of radius 200 AU with Keplerian rota- 
tion (iLaunhardt & Sarg ent 2001), surrounding a very young (ob- 
scured) low-mass T Tauri star. It is associated with a small 
bipolar near-infrared (NIR) nebula bisected by a dark extinc- 
tion lane at the posi tion and orientation of the edge-on disk 
(IStecklum et al.ll2004l) . The source is surrounded by an optically 
thin asymmetric envelope with a well-ordered magnetic field 
directed along RA. ~ 25 ° (iHenning et al.1 l200ll) . Furthermore, 
an Herbig-Haro object (HH494) was identified by Ha and 
S[II] narrow -band imaging, 6.15 arcmin northwest of CB 26 at 
RA. = 145° (IStecklum etal.l 120041) . The HH object is thus per- 
fectly aligned with the symmetry axis of the disk and the bipolar 
nebula. 

Here, we present new interferometric molecular line data of 
CB 26 and report the discovery of a small, well-collimated bipo- 
lar molecular outflow. The outflow has a peculiar velocity struc- 
ture that suggests it is rotating with the same orientation as the 
disk. 

2. Observations and data reduction 

Observations of CB 26 were carried out with the IRAM 
PdBI in November (D configuration with 5 antennas) and 
December 2005 (C configuration with 6 antennas). The two 
array configurations provided baseline lengths in the range 
16- 175m. The phase center was at 0^2000 = 04^59'^50.74^ 
^^2000 = 52°04'43.80''. Two receivers were used simultaneously 
and tuned single side-band (SSB) to the HCO+ (1-0) line at 
89.188526GHz and the ^^CO (2-1) line at 230.537984GHz, re- 
spectively, adopting the systemic velocit}0 of CB 26, vlsr = 
5.5 km s"^ In this paper, we present only the ^^CO (2-1) data on 
the molecular outflow. A more detailed analysis of the disk struc- 
ture, that includes other molecular line data, will be presented in 
a forthcoming paper. 

The quasar 0355+508 served as phase calibrator to determine 
the time-dependent complex antenna gains. The correlator band- 
pass was calibrated on 3C454.3 and 3C273, and the absolute 
flux density scale was derived from observations of MWC 349. 
The flux calibration uncertainty is estimated to be <20%. The 
primary beam size at 1.3 mm was 22''. The SSB system tem- 
perature at 230 GHz was in the range 250-450 K. Total band- 
width and effective channel spacing for the ^^CO (2-1) observa- 
tions where 12MHz (16kms"^) and 32kHz (0.25 kms"^), re- 
spectively. Zero-spacing data were obtained in September 2006 
with the IRAM 30 m telescope at Pico Veleta in Spain. They 
were combined with the interferometric data to help the decon- 
volution. However, the central velocity channels between 4.75 
and 6.75 km/s, which are completely dominated by the extended 
envelope, could not be well-restored and were therefore masked 
out for the subsequent analysis (Fig.O. 



The data were calibrated and imaged using the GILDASEl 
software. Channel maps were obtained with natural uv- 
weighting, resulting in an effective angular resolution of 1.47'' 
(synthesized beam FWHM). The r.m.s. in the channel maps is 
20 mJy/beam. 

New NIR K'-band images at 0'.'6 seeing and with 0'.'2 pixel 
scal e were obtained in J anuary 2008 with the Omega Cass cam- 
era dLenzen et al.l 11998 ) at the 3.5 m telescope of the Centro 
Astronomico Hispano Aleman (CAHA) at Calar Alto in Spain. 
These data are used here as complementary information to il- 
lustrate the morphology of the bipolar reflection nebula with re- 
spect to the edge-on circumstellar disk (Fig.[T]). Additional SMA 
interferometric data at 270 GHz to which we refer here (see also 
Fig.[T]) will be presented in a forthcoming paper (Launhardt et al. 
in prep.). 



3, The i^co (2-1) maps 

The ^^CO (2-1) channel maps and integrated intensity maps of 
CB 26 reveal a small, well-collimated, slightly asymmetric bipo- 
lar molecular outflow oriented perpendicular to the plane of the 
disk at P. A. = 147 ± 3°. This is illustrated in Fig.O which shows 
an integrated spectrum together with a set of binned channel 
maps. A complete set of channel maps is shown in Fig. ??, wich 
is available in electronic form at the CDS and can be obtained 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cds web.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 

The NW lobe of the outflow is pointing towards us (blue), 
the SE lobe is pointing away. This is consistent with the Doppler 
shift of the Herbig-Haro object HH494, located 6.15 arcmin 
northwest of CB 26 o n the axis of the molecular outflow (Fig.[Tl 
ISteck lum et al .1 120041) . At the central position, emission from the 
disk and outflow are superimposed. The total projected length 
of the bipolar molecular outflow is 14", which corresponds 
to ?^2000AU at 140 pc. However, the separation between the 
CB 26 disk and HH 494 of ^ 0.25 pc suggests the total extent of 
the outflow is much larger. 

Both the moment map of the ^^CO (2-1) emission (Fig. [3]) 
and the position- velocity diagrams (Fig. [4]) indicate the presence 
of a systematic velocity gradient perpendicular to the flow axis 
along the entire outflow. The orientation of this lateral velocity 
gradient is the same as that in the disk. The central position- 
velocity diagram across the disk shows the classical Kep lerian 
velocity pattern (Fig.Sl cf. ILaunhardt & SargentI l200ll) . The 
position- velocity diagrams at +3" and -4" away from the disk 
show clearly that the CO emission from the south-west side of 
the outflow axis (negative Ax) is red- shifted with respect to that 
from the north-east side (positive Ax) by about 2 kms"^ The 
amplitude of this lateral velocity gradient is comparable to that 
of the longitudinal velocity shift due to the bipolarity of the out- 
flow. Although this is probably a pure coincidence due to the 
high inclination angle with respect to the line of sight (see ? I4.1I) , 
it has certainly helped to disclose these kinematic features in the 
first place (see also discussion in ? 16.21) . The observed radial ve- 
locity field of CB 26 is very suggestive of a rotating outflow. To 
test this hypothesis and to derive more quantitative results, we 
develop in the following section an empirical model of the disk- 
outflow system and compare it to the full ^^C0(2-l) data cube. 



^ All radial velocity values in this paper refer to the Local Standard 

ofRest(LSR). ^ http://www.iram.fr/IRAMFR/GlLDAS 
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4. Modeling disk and outflow 

In order to reproduce the observed ^^CO (2-1) maps of CB 26, 
models of both the disk and the outflow are required. In this 
paper, we use a fixed disk model that was derived from dust 
continuum observations and develop an empirical model for the 
outflow. We then compute synthetic CO maps for the combined 
disk+outflow system and compare them to the observations (on 
the T^ij scale). Note, that in the following we use two different 
notations for the radius vector: "7?" refers to cylindrical coordi- 
nates, while "r" refers to spherical coordinates (see Fig. [6]). 

4.1. Disk model 

For the disk, we use the best-fit empirical model derived from 
multi-wavelength continuum observations of CB 26 (Sauter et 
al. in prep.). The basis of this model is a Keplerian disk with the 
following density distribution: 

\21 



Pdisk 



exp 



h{R) 



with h{R) = /zo 



(1) 



Here, R and z are the cylindrical coordinates and Rq = 100 AU. 
Furthermore, the model includes an optically thin extended en- 
velope to account for the scattered light observed at NIR wave- 
lengths. The disk is assumed to be passive and is heated by 
the central star only. Values of the free parameters (po, ap, ho, 
and ah) were derived by radia tive transfer simulations with the 
MC3D code (IWolf et al.ll2Q0l to fit the observed SED and the 
spatially resolved maps at JHK NIR-bands as well as at three 
millimeter wavelengths. In particular, the following parameters 
were derived: ap = 2.2, ho = 10 AU, and ah = 1.4. The 
disk mass remains somewhat uncertain, since the grain prop- 
erties, and hence the dust opacities, are not well-constrained. 
Assuming that this young disk (see ? I5.1I) has a "standard" in- 
terstellar gas-to-dust ratio, its total mass would be > 0.1 M©. 
This mass is consistent with the simple estimate derived di- 
rectly from the 1.3 mm contin uum flux of 190 mJy and us ing 
A^i.Smm = l.Ocm^g"^ of dust (lOssenkopf & Henninglll994h . a 
"standard" gas-to-dust ratio of 100, and <rd> ~ 20 K. The exact 
value of the disk mass has very little impact on our result, as the 
CO emission is optically thick. 

The disk was found to have an inner hole of radius ^ 45 AU. 
Note, that the 1.1 mm continuum image presented in Fig. [T] 
shows two peaks, which were already strongly suggesting the 
presence of an inner gap. The disk modeling also indicates an 
outer radius of 200 AU. 

The adopted disk model was combined with line radia- 
tive transfer (LRT) analysis of the HCO+ (1-0) and i2CO(2-l) 
Keplerian velocity profiles to derive the dynamical mass of the 
central star(s) and the inclination of the disk, = 0.5 ±0.1 M© 
and / = 85 + 4°, respectively. This mass is somewhat higher, 
but still comparable to the value of ^ 0.3 M© derived by 
iLaunhardt & Sargent (2001) from ^^CO (1-0) only. 

The ^^CO (2-1) emission of this disk was then computed and 
added to the outflow model described below, to be compared to 
the observations. We want to note that the parameters derived for 
the outflow do not critically depend on the details of the adopted 
disk model. A much simpler disk model, which we used first, 
lead to the same conclusions for the outflow. 



4.2. Outflow model 

As there exists no unified physical paradigm for molecular out- 
flows, their prescription still relies on empirical parameterized 



models adapted to the particular source. To model the CB 26 out- 
flow, we use a conical geometry of maximal size ro (see Fig. [6]). 
To reproduce the shape of the position- velocity diagram along 
the jet axis (Fig.[5]), we assume that the radial outflow velocity is 
a linear function of the distance from the star: V(r) = Vq • r/ro. 
Such "Hubble" law-type veloci ty fields have bee n observed in 
many molecular outflows (e.g. JLada & Fich 1996), but their ori- 
gin is not yet well-understood (see also Downes & Cabrit 2003b. 
We further consider a steady- state model: at any radial distance 
r from the central star, r^n(r)V(r) = r^noVo, where Vq and no are 
the radial outflow velocity and density at distance tq. The radial 
H2 density distribution is then given by n(r) = no • (r/ro)"^. For 
the distribution of the kinetic gas temperature in the outflow we 
assume a power law: T(r) = To(r/ro)~^. 

We assume that the outflow is rotating around its polar axis. 
In principle, one can consider different prescriptions for the ra- 
dial dependence of the rotation velocity of the outflow. We chose 
a simple prescription, Vrot ^ and modeled the following 
cases: ji = I (solid body rotation),// = -0.5 (Keplerian rotation), 
and yu = -1 (angular momentum conservation). However, the 
data did not allow us to reliably distinguish between these cases 
and to further constrain the parameter ji. Therfore, we adopted 
the physically most plausible scenario in this range. We assume 
that the rotation velocity is given by VrotW = Vkepi^^i) ' RlIR, 
where R is the distance from the outflow axis and Vkepi^L) is the 
Keplerian velocity at Rl (see Fig. [6] and discussion in ? l5.3b . This 
prescription is equivalent to the conservation of the angular mo- 
mentum along the flow, if all material was ejected from the disk 
at distance Rl with the tangential component of ejection velocity 
equal to VkepiRi)- 

Finally, we added a microturbulent velocity of 
O.lkm/s, which is a typical value for protoplanetary disks 
( Guillote au & DutrevI 11998) and assumed that CO is not de- 
pleted and has a constant abundance of X(CO) = 7.5 • 10"^ (see 
? l5.3b . Synthetic ^^CO (2-1) channel maps were p roduced using 
the LRT code URAN(IA) (iPavlvuchenkov et al.l 12007). taking 
as input both the outflow and the disk models predictions. The 
free parameters Vq, 7?^, ^0, ^0, and q of the outflow model were 
then derived by fitting the synthetic maps to the observations. 
The computed maps were convolved with a Gaussian beam of 
1.5'' and compared with the observed channel maps using the 
criterion: 



X 



pos,V 



rinod 



jobs 



pos,V pos,V) 



(2) 



Here, r^f y and rJ'J^ y are the modeled and observed intensities 

for a given position (pixel) and velocity channel, SF^^ is the 
noise level of the observed intensity in the channel maps, and 
summation is made over all positions and velocity channels in 
the observed map (N = Npos x Ny). Because of the large number 
of free parameters, we performed the -minimization in several 
2D cuts of the parameter space. This does not guarantee the ex- 
act localization of the minimum, but gives a reasonable approx- 
imation to it. This approach is justified since our goal was not 
to derive precise outflow parameters but rather to illustrate that 
the observed CO emission can indeed be described by a rotating 
outflow. 
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5. Results 

5. 1. The disk and central star(s) 

Although we have no direct indication for a double star or double 
jet in CB 26, the large inner hole in the disk of Rin ~ 45 AU (see 
S I4.1I) suggests the presence of a binary system. Assuming that 
this hole is due to tidal truncation by a binary system with ap- 
proximately equal mass es, the components should have a sepa- 
ration of a ~ 20 AU (see lArtymowicz & Lub ow 1994). As noted 
above, the Keplerian profiles of the line emission indicate a total 

stellar mass of 0.5 ± 0.1 Mp. 

In a recent paper, iGuilloteau et al.l (l2008l) have derived very 
similar properties (stellar mass, size of the central gap) for the 
HH 30 disk. Combining with the results of the optical jet study 
of IXnglada et al.l (l2007l) . they also infer the presence of a binary 
system with separation a ~ 15 AU, orbital period P ~ 80yrs, 
and total mass ~ 0.5 M©, although the mass ratio remains poorly 
constrained. 

Integrating the spectral energy distri bution (SEP) of CB 2 6 
yields a total luminosity of Lboi > 0.5 L© (IStecklum et al.ll2004) . 
This can only be a lower limit to the intrinsic luminosity since 
the edge-on view and bipolar geometry cause an anisotropic ra- 
diation field where only a fraction of the NIR photons are scat- 
tered into the line of sight. Higher intrinsic stellar luminosi- 
ties of up to ~ 2 L© are consistent with the model (but imply 
slightly different morphologies of the disk and envelope). Our 
model corresponds to Lboi = 0.9 L©. Nevertheless, this lumi- 
nosity range does allow an estimate o f the age of the s ystem. 
Using solar-metalicity PMS tracks of iSiess et al.l (l2000h and a 
star mass of 0.5 M©, we derive an age of 1 Myr, with lower and 
upper limits of 0.6 Myr and 2 Myr, respectively. Assuming two 
stars with 0.25 M© each, we derive a similar age range, but the 
lower limit shifts to 0.1 Myr. However, given the well-developed 
disk and the relatively thin remnant envelope, such a young age 
does not seem realistic. We can therefore safely constrain the age 
of CB 26 to be ^^1 Myr with an uncertainty of factor 2. 

5.2. The outflow 

Table[T] summarizes the parameters of the best-fit outflow model. 
The radial velocity field of the model can be compared to the 
observations (Fig. [3]): it reproduces well the mean radial velocity 
pattern (red and blue areas) and the S-like shape of the zero off- 
set velocity emission (green area). Furthermore, the computed 
position- velocity diagrams show the same radial velocity shifts 
across the outflow as the observed data (Fig.|4|. Although we can 
qualitatively reproduce all significant features which are indica- 
tive of outflow rotation, the model fit is not perfect. In particular, 
we can not reproduce the asymmetry of the outflow along its 
axis and the exact shape of the position- velocity diagrams. We 
believe that this is a drawback of our simplified and symmetric 
model. 

The mean H2 density in the outflow at ro = lOOOAU is 
2x10"^ cm"-^. The kinetic gas temperature at this radius is >15 K. 
The ^^CO (2-1) emission is moderately optically thick through- 
out most of the outflow. The derived radial outflow velocity at 
1000 AU of 10 km s"^ is comparable to that of HH 30 ( Pety et al. 
l2006), but has a large uncertainty (up to an order of magnitude) 
because of the uncertain inclination (see below). The rotation 
velocity is ~1 km s"^ at a distance of 100 AU from the ejection 
axis. 

The dynamical timescale of the molecular outflow (ratio of 
projected length and mean radial velocity, independent of the 



inclination) is Tdyn ~ 500 yr. This is much shorter than the esti- 
mated age of the source of ^ 1 Myr (S I5.1I) . However, the pres- 
ence of the Herbig-Haro object HH 494 at a much larger distance 
(Fig.[T]and Stecklum_etalJ|2004j) indicates that the total size, and 
hence the age, of the outflow is much larger than what we cur- 
rently see in the ^^CO (2-1) line. 

The total mass of the visible molecular outflow of ^ 
10"^ M© corresponds to about 1% of the disk mass. Although 
it is at the low-energy end, the CB 26 outflow fits well on the 
observationally established source luminosi ty - outflow momen - 
tum flux relation for protostars presented by lRicher et al.l (l2000h . 

5.3. Modeling uncertainties 

Interferometric channel maps of molecular lines from cir- 
cumstellar disks can be used to constrain the inclination 
angle and derive the kinematic mass of the central star 
(iGuilloteau & Dutrevl[l998l: ISimon et al.ll2000 i). However, there 
is a general degeneracy between disk inclination and star mass. 
While for face-on configurations the inclination angle can be 
well-constrained and the star mass remains uncertain, the situa- 
tion is opposite for edge-on configurations such as CB 26. Using 
geometrical arguments, we could however constrain the incli- 
nation angle of the CB 26 disk to be <89° (assuming the radial 
outflow velocity is < 100 km s"^) and >80° (from the aspect ratio 
of the 1.3 mm dust continuum emission). Therefore, we adopted 
a value of / = 85° ± 4°. In this inclination range, the central star 
mass is well-constrained by the Keplerian rotation field and we 
derive = 0.5 + 0.1 M©. However, the projection factor 1 /cos / 
has a total uncertainty of one order of magnitude, which leads to 
correspondingly large uncertainties in estimates of outflow ener- 
getics and timescales. 

Another degeneracy exists between hydrogen density and 
CO abundance. This would lead to corresponding uncertainties 
in the estimates of outflow mass and angular momentum. We 
cannot directly derive the CO abundance from our data, but we 
argue that the mean kinetic gas temperature in the outflow of 
> 15 K is high enough to rele ase all CO m olecules from the grain 
mantles into the gas phase (lB isschop~et a l. 2006) and assume 
maximum possible CO abundance. 

A more severe uncertainty comes from the fact that the small 
number of usable velocity channels (due to the effect of the en- 
velope; see §|2] and Figs.O IH and [5]) and the moderate SNR 
of the data did not allow us to well-constrain the radial depen- 
dence of the outflow rotation velocity, Vroti^^)- The assumption 
of Vrot ^ was made because it is the physically most plausi- 
ble profile within the range of profiles that are compatible with 
the data. However, we also want to stress that the other outflow 
parameters and conclusions we derive in this paper do not crit- 
ically depend on the exact radial profile of the outflow rotation 
velocity. 

6. Discussion 

6. 1. Alternative explanations for the velocity field of the 
outflow 

The previous sections presented the model of a rotating outflow, 
that was successfully fitted to the data. However, we also consid- 
ered a number of other configurations that could possibly lead to 
the observed peculiar velocity field of the CB 26 outflow: 

1 . molecular clouds which intersect with the line of sight and 
mimic a rotation-like velocity pattern. 
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2. two mis-aligned (non-rotating) outflows that are not spatially 
resolved, and 

3 . a precessing (non-rotating) jet / outflow. 

Given the relatively isolated location of the small globule 
CB 26 (see Fig.[T]), the clear jet-like bipolar morphology of the 
CO emission, and the single-component line shapes (see Fig.O, 
scenario (1) would require a very special spatial and kinematic 
configuration of the environment and seems extremely unlikely. 
The alignment of HH494 (P.A. = 145°) with the axis of the CO 
outflow (P.A. = 147 ± 3°) also argues against this scenario (see 
Fig.ffl). 

We tried to model two (spatially unresolved) outflows that 
are slightly misaligned with respect to the line of sight (scenario 
2). The suspected binarity of CB 26 could make this scenario 
appealing. However, until now, only very few close binary sys- 
tems are known to drive two jets; e.g., there is evidence that 
the quadrupolar outflow associated with HH 1 1 1 is driven by a 
close binary (separation < 40 AU in projection: [Rodriguez et al.l 
|2008l). Two outflows from the same central core (the so-called 
quadrupolar outflows) are usually observed only from sources 
with much larger separ ations ( e.g., a few th ousand AU in HH 288 
or CG30; G uethet aLll2QQH IChen et al.l 12008). Still, we have 
modeled this configuration and found that it can also produce an 
effective velocity gradient across the map of the two combined 
flows. However, reproducing the S- shape of the zero- velocity 
zone in the moment map (green area in Fig.O requires fine- 
tuning of both the relative orientation and the strength the two 
flows. The most significant differences to a rotating configura- 
tion occur near the driving sources and at low radial velocities. 
Unfortunately, both the dominating disk emission at the center 
and the missing low-velocity channels due to masking by the 
non-recoverable envelope emission prevent us from ultimately 
confirming or disproving this scenario as possible origin of the 
peculiar velocity field of CB 26. Note, however, that the agree- 
ment between the velocity gradient directions of disk and out- 
flow would be a pure coincidence in this scenario. 

Entrainment by an unseen precessing jet (scenario 3) could 
possibly also mimic a similar radial velocity field when the out- 
flow traces about 0.25 to 0.5 of the precession period. In this 
case, we would not observe rotation, but see gas that is acceler- 
ated along different directions at different times. However, lon- 
gitudinal gradients would dominate the radial velocity field over 
lateral ones. Since simulating this scenario would require a non- 
stationary model, we have only considered it as a thought ex- 
periment. Assuming that the precession of the jet axis is driven 
by tidal interactions between the circumstellar disk from which 
the jet is launched and a companion star on a non-coplanar or- 
bit, the shortest possible precession period in CB 26 would be of 
order 2000 yr (Terqueml l 19981: Anglada et al.j|2007l see their eq. 
14). Here we assumed equal companion masses, small inclina- 
tion angle, maximum ratio between circumstellar disk radius and 
c ompa nion separation of 0.3, and an orbital period of ^lOOyr 
(S I5.1I) . Comparing this to the dynamical timescale of the molec- 
ular outflow of 500 yr (? I5.2I) leads to the conclusion that at least 
the timescales are such that this scenario cannot be excluded to 
be at work for CB 26. It would, however, be worthwhile to test 
the precession scenario with an appropriate model. 

6.2. A rotating outflow - what distinguishes CB26 from 
other sources? 

Assuming the velocity pattern of the CB 26 outflow does result 
from a rotation of the outflowing material about its ejection axis. 



the data presented in this paper would be the first observational 
signature of this phenomenon. 

If the angular momentum transfer from the disk by an (un- 
seen) rotating jet drives the rotation of the molecular outflow, 
the question that arises is why such a signature is so prominent in 
CB 26 while it has not been observed in other sources? Detecting 
the projected velocity gradients requires an nearly edge-on con- 
figuration, which is the case of CB 26 but also of other sources 
that have been studied. Some specificities of CB 26 must there- 
fore be invoked, but these are difficult to spot without being too 
speculative. 

T he comparison between the CB 26 and HH 30 (iPety et al.l 

'20061: iGuilloteau et al.ll2008 ») cases illustrates that difficulty. By 
many aspects, these two sources are very similar: both systems 
are seen almost edge-on, with inclination larger than -84°; the 
disks are in Keplerian rotation around a -0.5 M© system and 
have an ~40 AU inner gap; and both sources are driving compact 
molecular outflows. Binarity seems the most simple explanation 
for bot h the mm and optical appearance of HH 30 ( Anglada et al] 
2007i: IGuilloteau et al.l ii2008). while the binarity status of CB 26 
is still speculative. If CB 26 is indeed a binary system, it is likely 
to have similar characteristics to those of HH 30, as the masses 
of the central star(s) and the sizes of the cavities are similar. As 
for the CO outflow, the HH 30 flow has a somewhat larger open- 
ing angle (60°) than CB 26 (40°, see Table [TJ and presents only 
one lobe (which remains a puzzling feature, probably related to 
the environment of the protostar rather than to an intrinsic prop- 
erty of the ejection mechanism). More importantly, and despite 
a detailed analysis of the CO e mission, the HH 30 outflow does 
not show any sign of rotation (iPety et al .1120061) . If the velocity 
pattern observed in CB 26 is indeed rotation, this raises the ques- 
tion of why rotation is not observed in a quite similar object like 
HH 30? The most significant difference between the two sources 
is their age: CB 26 is likely to be younger, since it has more 
prominent envelope. However, there are no clear prescriptions 
or predictions on the evolution of jet/outflow rotation with time. 

More generally, we note that the angular momentum trans- 
fer from the central sources to the outflowing gas may be quite 
complex in a binary system - as suspected in CB 26. One could, 
e.g., speculate on the presence of two distinct parallel jets, that 
are creating one single outflow, but "stir" the entrained gas 
and disperse orbital ang ular momentum from the system (cf. 
iMasciadri & Ragall2002h . 

6.3. The relation between disk and outflow 

If we make the two simplifying assumptions that the outflow 
mass flux is constant over time and that the outflow is fed directly 
by the disk, it would take ^^0.5 Myr to dissipate the disk (see 
TableO. This compares well to the age of the system and, within 
the uncertainty range, to the ty pical dispersal time of circum- 
stellar disks around young stars (iHaisch et al.ll200ll: ICieza et aP 
.2007,) . We also note that this mass outflow rate is comparable 
to typical accretion rates for young T Tauri stars (iHartigan et al.l 
119951: iGullbring eran il998h . Assuming that the outflow indeed 
rotates, the average specific angular momentum of the outflow 
gas is somewhat lower or comparable to that of the disk mate- 
rial. If we make the same assumptions for the angular momen- 
tum dispersal as above, we derive a timescale of 1 Myr. The 
good agreement between both balance and timescale estimates 
supports the hypothesis that the CB 26 outflow is indeed rotat- 
ing and dispersing angular momentum from the system. In any 
case, it also shows that in CB 26 the outflow plays a significant 
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role in dispersing the disk, and thus terminating the conditions 
for planet formation. 

The rotation of the outflowing gas can also be used to put 
constraints on the size of the ejecting region in the disk, assum- 
ing we are observing gas that has been directly ejected from 
the disk, and not envelope material that has been accelerated 
by an underlying (unseen) jet. Using the simple assumptions 
that the outflow material is launched from a single point on the 
Keplerian disk and that angular momentum is conserved, i.e., 
R • Vrot(R) = Rl ' Vkep(^L), wc dcrivc from our;^^-fitting routine 
/?L ~ 30 AU. This is consistent with the simple estimation 



Rl = 



(3) 



which yields 7?l ~ 27 AU when we assume = 0.5 M©, R = 
100 AU, and Vrot = 1.1 km s"^ (see Table[I]). 

One can also use mag neto-centrifugal wind models (see re- 
view bv IShu et al.1 [2000) to constrain the outflow launch re- 
gion. In this theory, the magnetic field plays an essential role 
in extracting the energy and angul ar momentum from the disk 
and redirecting them to the outflow. [Anderson et al.l ([20031) have 
shown that there is a useful invariant that is conserved along a 
magnetic field line and that can be used to determ ine the radius 
Rml at which the material was ejected. [Pety et al. ( 2006) have 
shown that this relation can be written as 



7?MT — R\2 



^kep^rot 

y2 



2/3 



(4) 



where R is the distance from the outflow axis to the considered 
volume element, V^ep = ^^GM^^R is the Keplerian velocity at 
this position, and Vrot and V are the measured rotation and radial 
outflow velocities at that position. How ever, the Hubble law-type 
radial outflow velocity field (see ? I4.2I) implies that the launching 
radius derived with this equation depends on the distance of the 
catching point from the star. For eq.[4]to b e applicable, the catch- 
ing point should far away from the star ([Anderson et al .1 12003). 
On the other hand, the Hubble law describes the observed ra- 
dial outflow velocity field of CB 26 reasonably well at small and 
intermediate distances from the star, while the outflow veloc- 
ity appears to be more constant at larger distances (see Fig. [5]). 
If we therefore use y = 5kms"^ at r = 500 AU (see Tabled] 
and §[121), and further assume M, = 0.5 M©, R = 100 AU, 
and Vrot = l.lkms"\ we obtain Rml ~ 32 AU. This value is 
remarkably close to the values obtained from the simple esti- 
mates above. Using V = lOkms"^ at r = 1000 AU leads to 
^ML ~ 13 AU, which can be considered a lower limit to Rml- 
However, as all these estimates are based on a number of un- 
proven assumptions, we do not use them to draw conclusions on 
the physical launch mechanism. 

7. Summary and Conclusions 

We presented new millimeter-interferometric observations of 
the young edge-on TTauristar - disk system in the isolated 
Bok globule CB 26. Our ^2CO(2-l) line data at VfS spatial and 
0.25 km s~^ velocity resolution reveal a previously undiscovered 
collimated bipolar molecular outflow of total length ^ 2000 AU, 
escaping perpendicular to the plane of the well-resolved circum- 
stellar disk in Keplerian rotation. The disk has an inner hole of 
radius ^^45 AU, indicating the possible presence of a binary star 
system with separation ^^20 AU. The dynamic mass and age of 
the central obscured star(s) is 0.5 ±0.1 M© and ^ 1 Myr, respec- 
tively. 



The outflow shows very clear and consistent kinematic sig- 
natures that suggest it is rotating about its polar axis. We use an 
empirical parameterized outflow model combined with 2-D line 
radiative transfer calculations and -minimization in parameter 
space to derive a best-fit model and compute parameters of the 
outflow. We could show that the data are consistent with outflow 
rotation. This hypothesis is supported by the fact that disk and 
outflow seem co-rotating and fit together energetically. However, 
we could not ultimately exclude jet precession or two misaligned 
flows as possible origin of the observed peculiar velocity field. 

Assuming that the CB 26 outflow is rotating, we computed 
and compared mass, mass flux, angular momenta, and angular 
momentum flux of disk and outflow and derived disk dispersal 
timescales of 0.5 ... 1 Myr. Although our simplifying assump- 
tions and the model uncertainties do not allow us to draw more 
precise conclusions on the outflow formation and disk disper- 
sal mechanisms, this result confirms the typical disk dispersal 
timescale derived on statistical arguments. It also shows that the 
CB 26 outflow must play a significant role in dispersing the disk. 

Ultimately testing the proposed scenarios that could explain 
the peculiar radial velocity field of the CB 26 molecular outflow 
will most likely require to combine the radiative transfer calcu- 
lations with hydrodynamic models of jet formation and entrain- 
ment in multiple sources. 
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Table 1. Parameters of the outflow model for CB 26. 



Parameter 


Symbol 


fixed 


fitted 


Outflow ^17P 




1000 ATT 








40 deg 




Positional angle 


PA 


30 deg 


_ 


Inclination 


i 


85 deg 




Microturbulent velocity 


^turb 


O.lkms-i 




CO/H2 abundance: 


X(CO) 


7.5 X 10-5 




H2 density at ro 


no 




2- lO^cm-^ 


Radial outflow velocity at ro 


Vo 




lOkm/s 


Temperature at ro 


To 




15 K 


Temperature index 


q 




0.5 


Rotation velocity at R =100 AU 


Kor(lOO) 




l.lkm/s 


Table 2. Derived parameters of outflow and disk. 


Parameter 


Symbol 


value 




Total outflow mass 


Mco 


1 X 10-3 Mo 




Mass flux 


Mco 


2 X 10-^ Mo 




Total outflow momentum 


Pco 


1 X 10-2 Mo 


kms~i 


Momentum flux 


P CO 


2 X 10-6 Mo 


kms~i yr-i 


Total angular momentum 


'^CO 


1.5 X 10^ Mokm^s" 


Angular momentum flux 




3x 10^ Mokm^s-iyr-i 


Specific angular momentum 


jco 


1.5 X 10^0 km^s-i 


Disk mass * 


^disk 


1 X 10-1 Mo 




Disk angular momentum 




3x 10^ Mokm^s 


Disk specific angular momentum 


Jdisk 


3x lO^^km^ 


s-i 



* The disk mass is derived from the disk model of Sautter et al. (in prep.). 
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Fig. 1. Overview o f the CB 26 region . The optical true-color image (left) is based on wide-field Ha (blue), [SII] (green), and I- 
band (red) frames (iStecklum et al.l l20Q4). The globules CB 25 and CB 26 as well as Herbig-Haro object HH494 are marked in the 
image. The zoom panel (right) shows a NIR K-band image of the bipolar reflection nebula (grey-scale, QffS resolution), overlaid with 
contours of the SMA 1.1 mm dust continuum emission from the disk (red, contours at 20, 38, and 55 mJy/beam) and the integrated 
^^C0(2-l) emission (0.6 to 12.5 km/s) from the bipolar molecular outflow (green, contours at 0.5, 1,2, 1.9, 2.9, . . . Jy/beamkm/s). 
Beam sizes are shown in Figs. E] and O The reference position is ^2000 = 04^59^^50.74^ ^2000 = 52°04'43.80'^ 
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Fig. 2. ^^CO (2-1) spectrum and binned channel maps of CB 26. Top panel: average ^^CO (2-1) spectrum of the central ±2". The 
central velocity channels between 4.75 and 6.75 km/s were masked out because they are dominated by the extended envelope. The 
integration ranges for the four binned channel maps (bottom panels) are indicated by vertical dashed lines and marked A, B, C, and 
D. Bottom panels: binned ^^CO (2-1) channel maps of CB 26, rotated by 30°, with integration ranges indicated by dashed lines in 
the top panel. Contour levels start at 0.15 Jy/beam km/s. The white cross marks the center of the disk. The ^^CO synthesized beam 
size is indicated as grey ellipses in the lower left corners. The two outer maps show the high-velocity Keplerian wings of the inner 
disk. The two inner maps are dominated by the blue (upwards) and red (downwards) outflow lobes, respectively. 




5 -5 5 -5 



Ax [arcsec] Ax [arcsec] 

Fig. 3. ^^CO (2-1) integrated intensity maps (contours) and mean velocity field (V^ moment map, color) of CB 26, rotated by 30°. 
White contours show the 1.1 mm dust continuum emission from the disk as observed with the SMA (contour levels same as Fig.[T]). 
The ^^CO synthesized beam size is shown as large grey ellipse. The smaller and darker ellipse shows the 1.1 mm continuum beam. 
Left panel: observations. Right panel: best-fit model for ^^CO (2-1). Dashed lines refer to the y-coordinate of the position-velocity 
diagrams shown in Fig.lH 
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Fig. 4. Position- velocity diagrams of ^^CO (2-1) perpendicular to the jet axis. The coordinate along the jet axis (Ay) is written in 
top right corner of each panel. Contours start at 60mJy/beam (3 cr). Left panel: observations. Right panel: best-fit model 
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Fig. 5. Position- velocity diagrams of ^^CO (2-1) parallel to the jet axis at Ax = 0". Contours start at 60 mJy/beam (3 cr). Left panel: 
observations. Right panel: best-fit model 
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Outflow 



_ R, - Launch radius 

Star ' 

Fig. 6. Schematic sketch of the outflow model. dV is the selected 
element of the outflow, r and R are the radial and cylindrical 
coordinates of the selected element, respectively. The element 
is launched from the disk at radius Rl. tq is the length of the 
outflow (see S I4.2I) . Matter flow lines (schematic) are shown as 
solid lines with arrows, while pure geometric connecting lines 
are dashed. 



